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ABSTRACT 

Force balance considerations put a limit on the rate of AGN radiation momentum output, 
L/c, capable of driving galactic superwinds and reproducing the observed Mbh — o~ relation 
between black hole mass and spheroid velocity dispersion. We show that black holes cannot 
supply enough momentum in radiation to drive the gas out by pressure alone. Energy-driven 
winds give a Mbh — & scaling favoured by a recent analysis but also fall short energetically 
once cooling is incorporated. We propose that outflow-triggering of star formation by enhancing 
the intercloud medium turbulent pressure and squeezing clouds can supply the necessary boost, 
and suggest possible tests of this hypothesis. Our hypothesis simultaneously can account for the 
observed halo baryon fraction. 

Subject headings: cosmology-black hole physics-galaxies:elliptical and lenticular-galaxies: formation- 
galaxies: evolution-quasars: general 

Introduction 



There is a consensus that the powerful Ac- 
tive Galactic Nuclei (AGNs) play a crucial role 
in shaping the general properties of galaxies (e.g., 



Silk fe Reeslll998l:ICroton et alfeOOdlSiik fc Norman 
2009) and clusters of galaxi es Ce.g.. lDalla Vecchia et al 



j(e.fi., L_ 

2004; iNusser et al.l 120061: iMcNamara fc Nulsen 



o + 



20071 ). AGNs are powered by accretion onto super- 
massive black holes believed to reside at the cen- 
ters of most galaxies. An indication of the galaxy- 
black hole connection is the remarkable correla- 
tion between the black hole mass, Mbh, and the 
velocity dispersio n, a, of the spheroid al galactic 
components (e.g.. iGuftekin et ah 20091) . Any suc- 
cessful model for galaxy formation must provide 
an explanation of this correlation. Self-regulated 
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Fig. 1— Plot of r)M B Hc/(M g o-c) (for n = 0.1) 
versus a from the sample of Gultekin et al 09. The 
horizontal line is O.IMbh = M g a/c. Black hole 
momentum falls short of the required momentum 
for most galaxies. 
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black hole gr owth offers a natur al explanation for 
this relation ( Silk fc Reed 1998 ). Both radiation 
pressure and mechanical outflows deposit momen- 
tum into the protogalactic gas. If this results in 
a wind, fo r ce balance argumen t s fdFabianl 1999; 



Kin d l2003t iMurrav et alT l2005: I Thompson et al 
20051 ) but see (lSokerll2009» ) lead to the conclusion 



that winds driven by pressure of radiation from a 
central black hole can suppress the collapse of gas 
and hence regulate the growth of the black hole. 
However the available momentum is, we show, in- 
sufficient to give the required normalization of the 
Mbh — cr scaling (Section 2). 



Th e original self-regulation argument of lSilk fc Rees 
(|1998l ) relied on energy balance: AGN activity 
heats the galactic gas reservoir above the virial 
temperature, generating galactic winds and even- 
tually terminating gas accretion onto the black 
hole. However, energy-driven winds suffer strong 
radiative cooling losses: while the radiation heats 
the gas nearby the black hole, the gas expands 
but cools rapidly, making the process inefficient 
(Section 3). Our preferred solution is to introduce 
positive AGN feedback via triggered star forma- 
tion. We argue that this simultaneously resolves 
three problems: the required order-of-magnitudc 
boost in the Mbh — cr scaling (Section 4), the 
enhanced specific star formatio n rate in massive 
galaxi es (addressed elsewhere bv lKhochfar fc Silk 
(2010), and the shortfall in the halo baryon frac- 
tion (Section 5). 

2. Can radiation momentum-driven winds 
expel halo gas? 

A luminosity L/c = M g g(r) balances the gravi- 
tational force applied on gas mass M g by the dom- 
inant dark matter (DM) in a protogalactic halo. 
For isothermal DM and gas distributions we get 
g(r) = 2a 2 /r and M g (r) = f s r 2 g(r)/G where a 
is the velocity dispersion and / is the gas mass 
fraction. The force balance condition, with the 
simplifying assumption t hat the gas mass M a (r) 



lies entirely on the shell (jMurrav et al.ll2005l) . 



yields a minimal luminosity 



L = 



4/X 
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(1) 



If we assume that L is proportional to Mbh as 
for the Eddington luminosity, L Edd — AttcGMbhiti^/ct^ 



this condition translates to 
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very close to the observed Mb h — cr relation 
M B H, 8 /cr 2 4 - 24±0A1 = 10 012±0 08 dGiiltekin et al. 
2009), where M BH ,s = Mbh/1O 8 M and cr 2 = 
cr/^OOkms- 1 . 

These arguments do not take into account the 
lifetime of the AGN. The following considera- 
tion shows that black holes obeying the observed 
M bh — cr relation cannot generate enough energy 
in radiation in order to drive the gas out of the 
protogalactic potential well by radiation pressure. 
The work, L(r e — r)/c, done by radiation pres- 
sure in moving the gas from r to r e must be 
sufficient to bring the gas to the escape speed, 
Ve(fe) = y/— 2</>(r e ), at r e . Energy conservation 
then demands 



'<r e 



r) > ^M g v 2 e (r) , 



(3) 



where v e (r) = y/— 2<p(r) is the escape speed at 
r. The total energy radiated by the AGN during 
this process is Lt where r = f ° dr/v(r) < (r — 
r e )/v e (r). Therefore, the energy radiated by the 
AGN must be greater than 



E,, 



L(r e -r)/v e (r) = -M g v e (r) (4) 



For an isothermal sphere truncated at the virial 
radius we have V 2 (r) = 4cr 2 (l + ln -^) so that 
E m > Mgccr. The total energy rjM B HC 2 (r) ~ 0.1 is 
the efficiency factor) that could be extracted from 
the black hole must therefore satisfy, 



T]M BH C > M g a 



(5) 



as the condition for momentum-driven winds to be 
able to unbind galactic gas. The virial gas mass 
in a halo is M g = 7.3 x 10 11 cr 2 3 (/ (5 /0.1)/ft,(z) M Q , 
where h(z) = H /H(z) and we have used H = 
71kms _1 /Mpc. Therefore, 

tjMbhc > 5 x 10 8 cr 2 4 {fj0.1)/h(z) M Q c . (6) 

Using the observed Mbh~ct relation, we find that 
the black hole cannot unbind the gas by momen- 
tum even for a — 400km/s at z = 2. The shortfall 
is about an order of magnitude at z = 0. 
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Complementary arguments a re recently given 
by lAnderson fc Bregmanl ( 2010h . These authors 
confirm the ubiquitous baryon fraction deficiency 
in galaxies. They compare baryon fractions in 
galaxies with varying bulge-to-disk-mass, and 
demonstrate that the presence of a supermassive 
black hole does not result in a reduced baryon 
fraction. 

In figure Q] we plot the ratio 0.1MBHc/(M g a) 
versus a from the ga laxy sample in Table 5 of 



Gultekin et al 



(|2009h . As a proxy for M g , the 
total gas mass to be ejected, we conservatively use 
the V-band luminosities of galaxies in the sample 
multiplied by a stellar mass-to-light ratio of 4. For 
the majority of galaxies in this sample the max- 
imum momentum t]Mbhc falls short of M g a by 
a factor of a few: half the sample falls short of 
the required threshold by a factor of at least 4, a 
quarter by a factor of 10. The figure also implies 
that the ratio depends weakly on the velocity dis- 
persion. 

Multiple scatterings can modulate the momen- 
tum delivered to the shell by the radiation but fail 
to resolve the momentum shortfall. Here is our 
reasoning. In the scattering case, force balance 
implies that 



pg{r) = up 



L 



4ircr 2 

where p is the gas density. Assuming g(r) is con- 
stant over the shell, we integrating this relation 
over the gas shell to give gM g = tL/c, where r is 
proportional to a in the cosmological setting. As- 
suming L is proportional to Mbh and since gM g 
is proportional to <r 4 , this force balance relation 
implies that Mbh °c a 3 . Modulation by multiple 
scattering is therefore inconsistent with observa- 
tions because r scales with a. 

To obtain the correct relation Mbh oc er 4 , or 
possibly steeper, we consider a photon delivering 
its momentum in a single encounter with a gas 
particle. This is achieved if the gas surrounding 
the black hole is mostly neutral, which is plau- 
sible if the cooling time is short in the central 
regions. Predominance of molecular gas in the 
nuclear regions is inferred indirectly from dust, 
CO and star formation observations. The opti- 
cal depth for dust also scales with surface density 
(jThompson et alj|2005l) . 



3. The case against energy-driven outflows 

Energy-driven outflows hav e been argued to 
give the wrong scaling relation ([Silk fc Reesl ll998) 
for the observed Mbh — c relation, although a re- 
cent reanalysis favours th e originally predicte d 
Mbh ~ & 5 dependence ( Graham et al.l l2010h . 
However there is a more fundamental difficulty 
with energy-driven outflows. The gas initially 
cools by Compton sca ttering wit h the radiation 
emitted by the AGN (jKinsd 120031) . This cooling 
indeed is important in the central region but lasts 
only as long as the AGN is active. Radiative cool- 
ing plays an important role in the AGN-galaxy in- 
terplay over longer timescales. Radiative cooling 
is extremely efficient in small halos ( 5.10 11 M Q ), 
where the co oling radiu s can even exceed the 
virial radius (Silk 19771: iRees fc Qstrikerl Il977: 



White fc Reesl 119781 ). To explore the role of 
radiative cooling, we have simulated the feed- 
back effects in spherical systems using a one- 
dimensional Lagrangian hydrod ynamical code (see 
iNusser fc Pointecouteaul (|2006l ) for numerical de- 
tails.) In these simulations, the dark matter is 
assumed to reside in a static isothermal spheri- 
cal halo truncated at the virial radius, R v , and 
the gravity of the gas is ignored. The gas is rep- 
resented by 250 shells, which are equally spaced 
between r = .1R V and R v , and are initially in hy- 
drostatic equilibrium in the gravitational potential 
of the dark matter, with zero external pressure at 
r = R v . AGN feedback is introduced as heat 
added to the innermost shell over a time-scale of 
5 x 10 7 yr, which is shorter than the dynamical and 
the radiative cooling time-scales. The equilibrium 
radiative cooling cu rve with metallicity of one- 
third solar as given in lSutherland fc Dopita ( 1993f ) 
is adopted. Figure [2] shows the results from the 
simulations for two halos with a = 100 kms^ 1 and 
er = 300 kms -1 . The explosion energy is taken 
to represent the AGN feedback, and is equal to 
the gas potential energy in absolute value. The 
integrated AGN energy delivered to the system 
is set equal to the initial potential energy of the 
gas, about 1.4 x 59 ergs and 3.4 x 61 ergs for 
cr = 100 kms' 1 and 300 kms -1 , respectively. The 
results in the left panel (no cooling case) demon- 
strate that for a = 100 kms -1 without cooling, 
the gas is ejected from the system. For this smaller 
halo the radiated energy is so large that the total 
energy of the system at the final time is nega- 
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tive and the gas falls back onto the halo, as seen 
in the middle panel. Only for massive galaxies, 
cr = 300 kms -1 , right panel, does cooling play a 
lesser role, so that the feedback actually manages 
to unbind the gas out of the halo. For the effect 
of the cooling to be negligible, a self-regulating 
mechanism must operate so that the gas is rapidly 
ejected before cooling becomes important. But 
such a contrived mechanism will likely produce 
a complicated functional form for the Mbh — & 
relation which should reflect the dependence of 
the cooling curve on a and the details of how the 
energy is ejected. Our discussion of the relevance 
of coolin g demonstrates that mo mentum-driven 
outflows (|Fabian|[l99l lKing||2005l) dominate over 
all mass scales of interest and, as we have shown, 
possibly give the correct scaling. Such flows are 
inevitable, as radiative cooling dominates over a 
longer time scale, especially in smaller systems. 
However, they cannot account for the observed 
normalisation of the Mbh — cr relation. 

4. The role of star formation triggering by 
AGN momentum-driven winds 

If AGN and SNe fail, admittedly for different 
reasons, to drive the required outflows, then we 
now argue that the combination may provide an 
ideal solution. For a black hole emitting at the Ed- 
dington luminosity, force balance yields Mbh ~ 
cr 4 , close to the observed scaling (see eq. [5]) . The 
problem as we have seen is that the black hole can- 
not supply sufficient energy to suppress its own 
growth by expelling gas by radiation pressure, i.e. 
the black hole does not operate for long enough. 
The shortfall is a factor of a few. This conclu- 
sion i s also sustained by recent observational data 
(e.g.. iDunn et al.l l2010h . We propose here that a 
radiation momentum-driven wind triggers a star- 
burst which joins forces with radiation to drive 
gas out of the protogalaxy. In this way, radiation 
momentum-driven winds, under certain physical 
assumptions, could yield a near Mbh ~ f 4 scal- 
ing although they are not solely responsible for gas 
expulsion. 

If the triggered starburst is to aid in unbind- 
ing the gas by heating it above the halo's virial 
temperature ,then NsnEsn ~ 0.5M g a 2 , where 
Esn ~ 10 51 ergs is the mechanical energy release 
per SNII and N$n is the number of SNII in the 



starburst. Further, the whole process of gas re- 
moval has to be short enough so that radiative 
cooling is not important. Since M g ~ a 3 and if 
we assume N$n ~ Mbh, we get Mbh ~ c 5 . This 
corresponds to the observed scaling in the energy- 
driven wind case and allows an acceptable nor- 
malisation. Alternatively, to achieve a cr 4 scaling, 
our preferred model invokes considerations of mo- 
mentum balance. The starburst must supply the 
"missing momentum" of a factor of a few of the 
black hole's t/Mbhc- The momentum boost works 
as follows. A supernova remnant (SNR) conserves 
energy until the swept- up shell mass decelerates 
to a velocity below v s w 400km/s (jCioffi et al 



1988), below which momentum is approximately 
conserved (in a uniform medium). The result- 
ing inefficiency of SNR energy input is of order 
Vci/vsn, where v c \ is the interstellar cloud veloc- 
ity dispersion and vsn ~ E$n /[msNVs], where 
Esn ~ 1.5 x 10 51 ergs is the initial SN explosion 
energy and msN ~ 150M Q (for a standard IMF) 
is the mass in stars required to form a supernova 
of type II. This inefficiency, of order a percent, 
is confirmed for more realistic conditions by nu- 
merical simulations of S N-driven galactic winds 
(|Dubois fc Tevssier|[2008h . 

The momentum input is Esn/v s per super- 
nova, and momentum balance for N$n supernovae 
gives 

NsnEsn/v s = M g a . (7) 

Since the gas mass M g ^ a 3 this yields a boost 
satisfying the desired dependence on a. The mo- 
mentum boost NsnEsn fas must be larger than 
t]Mbhc by a factor of a few. To obtain this, we 
assume that the black hole growth rate is equal to 
the Eddington-limitcd accretion rate and that the 
star formation luminosity is equal to the Edding- 
ton luminosity. This is a re asonable appr oxima- 
tion at zeroth order: in fact Netzer (2009) shows 
that there is a small tilt in the logarithmic rela- 
tion for the star formation luminosity dependence 
on the luminous AGN that concern us here, with 

t oc r Q.g 

where the bolometric luminosity is the sum of 
star formation and AGN luminosity. Our ne- 
glect of this tilt leaves the boost factor indepen- 
dent of mass and preserves the cr 4 dependence 
of the Mbh,o correlation. Therefore, we sim- 
ply assume that the total star formation energy 
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Fig. 2. — The effect of radiative cooling on the ejection of gas by a central explosion in isothermal halos. 
Left, middle and right columns, respectively, correspond to a = 100 kms -1 without cooling, a — 100 kms -1 
with radiative cooling, and a = 300 kms -1 with cooling. Solid red lines represent the hydrostatic initial 
conditions. Blue dots represent the systems at some intermediate time. Black dots correspond to the final 
time of the runs at 1.3 Gyr after the explosion. Note that this 1-d simulation does not describe fragmentation 
of shells. 



is E„ = f*r]MBHC 2 , where /* is a factor of or- 
der unity. The stellar mass associated with E* is 
M* = E*/(e nuc c 2 ) = /,(??/0.1)10 10 M Bffi8 , where 
e nuc s=s 10 -3 — 10~ 4 is the thermonuclear burning 
energy efficiency for a massive star. The associ- 
ated mechanical energy produced by SNII in the 
starburst is 



f*E*E S N , 
— = L.«Wff c 



NsnEsn/vs 

(8) 

Also Nsn = Mz/msN, where itisn ~ 100M Q is 
the mass in stars formed per SNII and is weakly 
IMF-dependent. The boost factor, 



f*E, 



SN 



£nucm S NCV s 



(9) 



amounts to an order of magnitude for /* = 1, 
E SN = 10 51 ergs, e nuc = 3.10- 4 , m SN = 15OM 
(expected for e.g. a Chabrier IMF), and v s — 
400km s _1 . Of course the question remaining is 
why the boost factor should be only weakly de- 
pendent on black hole mass. This requires more 
than a universal mass function of gas clouds since 
the AGN outflow pr essure enters. But this might 
work, see eq. 85 in ISilk fc Norman (|2009f) : the 



AGN-driven supersonic turbulence velocity dis- 
persion is found to depend only logarithmically on 
AGN properties, as also does the porosity which 
controls turbulent pressure. 

We now discuss the implications of our feedback 
models for further applications, notably to satel- 
lite abundances, intermediate mass black holes, 
and the baryon fraction. 

5. Feedback, satellite abundances and 
baryonic fraction 

AGN-triggered preheating is a plausible mech- 
anism for reducing satellite abundances in galaxy 
groups or around m assive early- type galaxies. 
Koposov et al. I (l2009h demonstrated that SN feed- 
back plus reionization accounts for the luminos- 
ity function of the MWG and M31 dwarfs below 
10 8 M Q . However, there are more intermediate 
mass galaxies(namely LMC/SMC/M32/NGC205) 
in the MWG and M31 halos than are found 
in models tuned to fit the ultrafaint dwarf fre- 
quency. Models which fit the ultrafaint dwarfs are 
so efficient that they underproduce the massive 
dwarfs/intermediate mass galaxies. Conversely 
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models tuned to fit the massive dwarfs have ex- 
cessively inefficient feedback and overproduce the 
numbers of ultrafaint dwarfs. This problem seems 



Smith et all (2009) 



to be common to all SAMs. 
confirms that SAMs fail to resolve this problem 
for the abundance of intermediate mass galaxies. 
This data set confirms that massive galaxie s are 
overproduced in the models. Liu et al. ( 20091 ) also 
finds that all SAMs overpredict the number of 
satellites by at least a factor of two in the mass 
range 10 9 - 10 10 M Q . 

AGNs are commonly introduced at early epochs 
to account for the black hole-bulge correlation via 
the quasar feedback mode at early epochs. At late 
epochs, the AGN radio mode inhibits cooling of 
the dilute gas resulting from the earlier feedback 
process, keeps the gas hot, resolves the galaxy 
luminosity function bright end problem, and ac- 
counts for the red colours of massive early type 
galaxies. We point out here that AGN feedback 
in the radio-quiet mode may also account for the 
suppression in numbers of dwarf satellite galax- 
ies. Feedback from AGN in the host galaxies pre- 
heats the halo gas that otherwise would be cap- 
tured by satellites. However, suppression of the 
formation of intermediate mass satellites of the 
MW, and more generally, late-type galaxies with 
small bulges may not be efficient because of the 
low masses of their central black holes. We have 
no solution to this problem, more generally as- 
sociated with the large observed frequency of es- 
sentially bulge-less thin disks, other than to sug- 
gest patchy accretion of cold gas must play a role 
in thin disk formation at late epochs (but see 
Peebles fc Nusserl i201(t \). Once the potential well 
of a massive galaxy has developed, SN do not eject 
gas, although they may drive interstellar turbu- 
lence and fountains. Gravitational heating does 
not work en route to forming the potential. All 
that is left is vigorous activity in the MW assem- 
bly stage. This phase may plausibly involve feed- 
back from intermediate mass black holes, which 
are believed, at least by some, to be ubiquitous. 

The hypothesis that intermediate mass black 
holes are formed generically during the hierarchi- 
cal build-up of galaxies may possibly provide a 
radical solution to the baryon fraction problem 
via the momentum-driven outflows that we are in- 
voking. Theoretical arguments suggest that one 
pathway towards building up the central SMBH 



is via mergers of intermediate mass black holes 
(IMBH) during the hierarchical merging evolution 
of dark matter halos. It is assumed that substruc- 
tures develop IMBH at early epochs, contempo- 
raneously with first star formation. Simultane- 
ously, another major problem is resolved, that of 
the baryon fraction, via preheating or ejection. 
This is seen to be low in low mass and in mas- 
sive galaxies (McGaugh et all 2009). If satellites 
form in a secondary manner, preheating reduces 
their baryon fraction. If satellites formed first, the 
quasar mode will sweep the gas out of the galaxy. 
This is achieved by a combination of AGN out- 
flow momentum plus induced SN feedback. The 
gas subsequently stays out via the quiet mode of 
AGN feedback from the central host and other ac- 
tive galaxies. 

Within the standard paradigm, the observed 
baryon fraction of the MWG with its small black 
hole is only explained if there are sufficient IMBH 
in the satel lites to drive out the b aryons at early 
epochs (see Keselman et al.l ( 2010() for an alterna- 
tive non-standard explanations). For galaxies with 
massive spheroids, the central BH strips the satel- 
lites, reducing the number of IMBH, but provides 
enough feedback to eject the baryons. 

Baryons must be ejected in order to begin with 
the primordial baryon abundance, from galaxies 
of all mass scales, as well as galaxy groups. Only 
at clusters scales is the baryon fraction conver- 
gent. For dwarfs this is not a problem, but for 
typical galaxies, such as the MW, SN feedback 
cannot be responsible for baryon loss. Rather the 
baryons are recycled via the halo. AGNs provide 
the only energy source capable of accounting for 
baryon ejection. 

6. Discussion 

The primary aim of our paper is to highlight the 
scaling relation normalisation problem for AGN 
feedback, and to propose a possible solution in- 
volving AGN-triggered star formation. Positive 
feedback may have important ramifications for 
star formation at high redshift, and is inevitably 
followed by gas outflows driven by both AGN and 
supernovae, along with concomitant quenching of 
star formation. 

From the data plotted in figure 1, the momen- 
tum boosting by the starburst is a factor of a few. 
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This will naturally yield the dispersion in the re- 
lations given the nature of the boost, e.g. by BH 
outflow triggering of SNII. The points that lie low 
in the momentum condition had a larger boost, 
and this would lead to a prediction that the resid- 
uals in Mbhc vs Mg<j should anticorrelate with 
SNII tracers in chemical evolution, e.g. the bulge 
a/Fe. 

Small galaxies which formed stars before host 
galaxy AGN onset will survive. They should be 
seen as a bump in the galaxy luminosity func- 
tion ( GLF), analogously to what is seen in the 
MW (IKodosov et all l2009h and in the K-band 
GLF (jSmith et al.ll2009h . These galaxies are dis- 
tinguishable by being older and more metal-poor 
than their AGN-modulated successors which are 
primarily either low mass satellites or massive 
ear ly- type galaxies . For the MW, the failure of 
the [Koposov et aL I (l2009h model tuned to the nu- 
merous ultrafaint dwarfs to account for the ad- 
mittedly sparse numbers of massive dwarfs is con- 
sistent with the lack of a large BH in our AGN 
feedback model. Feedback from IMBH can re- 
solve this problem. We have suggested that the 
IMBH in wCen may be an example of a popula- 
tion of halo IMBH that could have provided the 
additional feedback needed to both allow the LMC 
and similar dwarfs to form and not simultaneously 
overproduce the faint dwarfs. Such IMBH could 
easily, during an active accretion phase, have pro- 
duced enough momentum to have swept the resid- 
ual gas out of the outer halo. 

Globular clusters are plausibly the most visi- 
ble surviving component of the first generation 
of substructure. That they might have a direct 
connection to IMBH is weakly supported by the 
possibility that one of the most massive globular 
clusters, cjCen, might contain an IMBH. Another 
hint of a connection with globular clusters may be 
present in the apparent correlation between black 
hole mass and mass of the host galax y globular 
cluster system (jSpitler fc Forbes! 120091 ) . A varia- 
tion on this relation has recently been found that 
relates blac k hole mass to the number of globu- 
lar clusters ( Burkert fc Tremainell201oj ). Numeri- 
cal simulations find that the SMBH-cr scaling re- 
lation c an be preserved by hier archical mergers of 
IMBH (| Johansson et al.ll2009l) . This lends sup- 
port to the possibility that globular clusters may 
serve as a proxy both for IMBH and for dwarf 



galaxies, and therefore provide a possible witness 
to the required baryonic cleansing role of satellites 
by IMBH in our model. 

We thank Noam Soker and Avi Loeb for valu- 
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